On Critical Exponents and the Renormalization of the Couphng Constant 
in Growth Models with Surface Diffusion 
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It is shown by the method of renormahzed field theory that in contrast to a statement based on a 
mathematically ill-defined invariance transformation and found in most of the recent publications on 
growth models with surface diffusion, the coupling constant of these models renormalizes nontrivially. 
This implies that the widely accepted supposedly exact scaling exponents are to be corrected. A 
two-loop calculation shows that the corrections are small and these exponents seem to be very good 
approximations. 
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There has been great interest in nonequilibrium growth 
processes of surfaces and interfaces in the last ten years 
(for general overviews on the work in this area of sta- 
tistical physics and the related literature see the articles 
of Halpin-Healy and Zhang , Krug |^,|| and the book 
of Barabasy and Stanley Q). Recently deposition pro- 
cesses in which surface diffusion constitutes the dominant 
relaxation mechanism have been intensively studied. Be- 
sides the simulational work on many different models, 
stochastic differential equations for the surface growth 
with surface diffusion arc proposed and studied by scaling 
and one-loop momentum-shell renormalization methods 
in order to find out the universal scaling behavior. An 
apparently mathematically ill-defined invariance trans- 
formation leading to supposedly exact scaling exponents 
plays a crucial role in this analytical work. This transfor- 
mation and its consequences have been taken over, appar- 
ently without careful examination, in most of the related 
publications and overviews , and are now widely ac- 
cepted in the community of surface-growth physicists. 

In this letter I explicitly show that the widely accepted 
invariance is indeed wrong and the supposedly exact ex- 
ponents are to be corrected. I will use the systematic 
method of renormalized field theory (see e.g. [^|-||) and 
calculate the corrections in an e-expansion to second or- 
der. 

I consider a system evolving by diffusion of atoms at 
a surface of (suitably scaled) height s (x, t) over a d- 
dimensional hyper-plane. Thus the deterministic part 
of the equation of motion for the stochastic dynamical 
variable s has to be of the conserved form 
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TS + V^S+^ {Vsf 



c 



(C (x, t) C (x', t')) = 2A-1 (-V^)" S (x - x') S {t - t') . (1) 

Here the parameter a = 1 in the fully conserved case orig- 
inally studied by Sun, Guo, and Grant |^] in their first 
exploration of non-KPZ universality classes (with r = 0). 
If the noise is not conserved (shot noise) the parameter 



a = 0. This case of (|^) was introduced in the first pa- 
pers on ideal molecular beam epitaxy (MBE) ||To|-[l^. In 
the Langevin equation (^) A is a kinetic coefficient and g 
is the coupling constant. The behavior near the critical 
point T = is of primary interest. For r > the equa- 
tions (0) are eventually driven by renormalization to the 
linear Edwards- Wilkinson equations both with noncon- 
served or conserved noise. 

For a generalization of the Galilean invariance (tilt 
symmetry) of the KPZ equation p3[ Sun, Guo, and 
Grant introduced as a symmetry of (111) the coordinate 
transformation 



x ^ x — XgtaW 



(2) 



which appears to be mathematically ill-defined because 
it mixes coordinates and linear differential operators in a 
coordinate substitution rule (the questionableness of the 
rule (^ shows up in a definition of the transformation of a 
product of functions /i (x) /2 (x) = /2 (x) /i (x)). Based 
on the transformation (|^) the authors state that the cou- 
pling constant g changes only trivially under renormal- 
ization and derive supposedly exact scaling exponents. 
The ill-defined transformation (||) has been used for vari- 
ous equations of the type (|l|) to derive critical exponents. 
Only few authors express some doubt concerning the use 

of d) 

In the following I show by an explicit two-loop cal- 
culation that the coupling constant does indeed change 
non-trivially under renormalization. Therefore the sup- 
posedly exact exponents are to be corrected by a quantity 
S: 



Z = 4: 



(3) 



where s = — d is the deviation from the upper crit- 
ical dimension d^. The roughening exponent x ^^'^ the 
dynamic exponent z are defined by the structure func- 
tion (|s(x,t) -s{0,t)f) ~ |x|^^/(V|xr), where / is a 
scaling function. I calculate the correction 6 to order 
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and find that it is always a small but nonzero positive 
quantity which has not been detected in simulations up 
to now. 

To develop a renormalized field theory I introduce 
the dynamic functional [^-f8^ corresponding to the 
Langevin equations ([^) 



J = j (fxdt\^\-^ls + TVT^s 



V^. + f (V.) 



(4) 



Here 's denotes the response field (the MSR auxil- 
iary variable). Now all response and correlation func- 
tions can be calculated by path integrals with the weight 
exp (— [s^, s]) . The expansion of the path integrals 
about the Gaussian part of J generates the (diagram- 
matic) perturbation theory which has to be regularized 
and renormalized because the series are UV-divergent at 
the upper critical dimensions that are given by dc = 4 
for a = and dc = 2 for a — 1. The propagator 
G'(q,t) = e{t) exp(-K(g)i) with k (g) = Xq^ {t + q^) 
and the correlator C(q, t) = X (?) cxp (— k (g) |<|) with 
x{l) = g^*-""^-* {t + q^^ are considered as functions 
of the momentum and time. The step function 9 (t) 
is defined with 6 (0) — 0. The vertex is given by 
^(qi;q2,q3) = q2 •q3w('7i) with v{q) = -\gq^. Here 
c\2 , qs are the ingoing momenta of the two s-legs and 
qi = q2 + ^3, denotes the outgoing momentum of the 
J-leg. Thus the diagrammatic expansions lead only to 
one-particle irreducible vertex functions with (at least) 
a factor at each (amputated) s-leg and a factor q at 
each s-leg. A look at the naive dimensions of these ver- 
tex functions shows that only primitive divergencies arise 
which correspond to the terms V^sV^s and V^s^(Vs)^ 
in the dynamic functional (^. Therefore the model is 
renormalizable by the scheme 



°s = Z^'^s, s^s^ Z-^l'^s, X^'X = ZX 
T ^ f = Z~ V, .g^ ^ 5^ = Aefi^Z^^ZuU. 



(5) 



is a suitable constant and fi is 
the usual inverse external length scale which makes the 
renormalized coupling constant u dimensionless. The two 
renormalization constants Z and Zu are determined by 
minimal renormalization. One easily gets (corresponding 
to the known 1-loop results) 



Z =1-2 



a-2 ' 



Oiu' 



(6) 



The incorrect assertion based on the invariance transfor- 
mation (^) is that Zu = 1 to all orders of u. In the 
following I give the derivation of the second order contri- 
bution. 

The calculation is greatly simplified due to the follow- 
ing properties of the diagrammatic contributions. Con- 
sider a diagram in the momentum-time picture. In 



this representation the restrictions which causality im- 
poses are more easily exploited than in the momentum- 
frequency picture; note in particular the ^^-function of 
the propagator. Due to causality the vertices are ordered 
from right to left with increasing time labels. Thus, a di- 
agram with n vertices consists of {n — 1) time segments. 
A Fourier transformation from time differences to fre- 
quencies then generates a dynamical factor 



^({q,^}) = 



l^UJa +^K(qfc) 



(7) 



for each time segment. Here uja is the sum over all 
external frequencies flowing into the segment from the 
right, and the set {qfc} are the internal momenta run- 
ning from right to left within the segment. In addition 
to the dynamical factors, each correlator contributes a 
factor X (q) and each vertex a factor V (qi; q2, qa)- Fi- 
nally one has to integrate over the internal momenta. 

I now consider one-particle irreducible diagrams with 
amputated legs as contributions to the vertex functions. 
For the determination of the UV-divergencies all exter- 
nal momenta and frequencies can be set to zero if the 
external momenta coming from the derivatives at the ex- 
ternal legs are included in the amputation procedure. To 
find the renormalization factor Z„, one needs the ver- 
tex function which results from the diagrams with am- 
putated legs corresponding to one V^s and two (Vs). 
I now discuss how certain diagrams cancel (this discus- 
sion is analogous to that in |20|). Consider the earliest 
vertex of the diagram. If this vertex has an amputated 
(Vs)-leg, its contributions coming from different permu- 
tations of the propagator (Fig. la) add to zero because 
q^ {—q) X (?) ~ q^^ (?) X (~9) — 0- it is essential for this 
cancellation that the vertex considered carries the ear- 
liest time argument, since otherwise the straightforward 
permutation of propagators is prohibited by causality. 
Especially, by this cancellation all one-loop vertex func- 
tions with zero external momenta vanish, and it follows 
that the coupling renormalizes only trivially to this order. 
If the earliest vertex is an internal one, the contributions 
coming from the propagator permutations are added up 
to the summed vertex (Fig. 2a) 



'S'(qi,q2,q3) = 



V"(qi;q2,q3) 

xiqi) 
^ 't^(q2;q3,qi) 

x{q2) 



^(q3;qi,q2) 

X(93) 



(8) 



with qi + q2 + q3 = 0. The factors result from the 
exchange of a propagator line with the corresponding cor- 
relator line. 

If in the diagrams under consideration the two vertices 
with the amputated (Vs)-legs are deleted, all possible di- 
agrams with internal earliest vertex (the only ones that 
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are nonzero) reduce to one kernel-diagram (Fig. 2b) with 
a contribution 



■f^(qi,q2,q3) 



Affgi^ V"(qi;q2,q3) 
2 2k (gi) 

x(gi)x(g2)x(g3) 

K{qi) + k((?2) + k(<73) 



S'(qi,q2,q3) (9) 



which is self-explaining (the first factor 1/2 results from 
the double connection of the two internal vertices). In- 
serting back the external vertices allows the discussion of 
a second cancellation mechanism. The insertions of an 
amputated vertex in a time segment with only two in- 
ternal lines add to zero if all external momenta are zero 
(Fig. lb) because qv (q) — qv {—q) = 0. Therefore only 
the insertions into the second time segment with three 
internal lines of the kernel-diagram (Fig. 2b) may lead to 
a nonvanishing contribution. The result of the addition 
of the insertions of one amputated external vertex with 
zero external momentum into a three-line segment (Fig. 
Ic) is the factor 



W(qi,q2,q3) 



qiv (qi) + q2« ((72) + q3'^ (g3) 
K{qi) + K{q2) + K{q3) 



(10) 



where qi -I- q2 + q3 = 0, and the denominator arises be- 
cause a new time segment is produced by an insertion. 
In the case of the KPZ equation the functions v (qi) are 
constants and it follows that these insertions add to zero: 
W = 0, a result which is easily generalized to insertions 
in all segments and has its origin in the Galilean invari- 
ance po| . But in the models considered here v (q) ^ q^ , 
and it is easily seen that W (qi,q2,q3) = only in the 
case of at least one vanishing momentum q^. In gen- 
eral W 7^ which eventually disproves all statements 
on the trivial renormalization of the coupling constant. 
The insertions of the two external vertices then lead to 
the tensor W (X) W. But after integration over the mo- 
menta the integral has to be a rotational invariant. Thus 
it must be proportional to the d-dimensional unit ten- 
sor and its prefactor can be determined by taking the 
trace. Therefore instead of the integral with the tensor 
W (g) W only the integral with the trace • W must 
be calculated to find the result of the insertion of the two 
external vertices into the kernel-diagram. Eventually as 
the contribution of all two loop diagrams to the ampu- 
tated vertex function with zero external momenta and 
frequencies I obtain the integrand 

2 

(qi,q2,q3) = ^-fs^(qi,q2,q3) W(qi,q2,q3)^ (ii) 

with qi + q2 + q3 = 0. The factor 2 arises from the 
permutation of the external amputated s-legs. 

Performing the integral over the internal loop momenta 

R ^ (27r)"^''/nrf''gi'^(qi +q2 + q3)-^(qi,q2,q3) one 

observes that the IR-regulator r > can be neglected 
in all terms with exception of the dynamical factor 



(k (qi) + K {q2) + K (q3))~ . Then the integration over 
the variable s = q^ + q2 + qi extracts the e-pole, and in 
the remaining part of the integral the dimension d can 
be safely set to dc- Up to trivial angle integrations the 
result R reduces to a finite two-dimensional integral with- 
out any singularities in the integrand. The elimination of 
the e-pole from the renormalized vertex function finally 
leads to a definitely nontrivial renormalization constant 



e 



Ba = 



0.0057406 
0.0167998 



if 
if 



a = 
a = 1 



(12) 



The RG fmictions are defined by 7 = (91nZ/i91n/i)Q, 
7ti = (91nZ„/01n/j,)p, and /3 = fj, (du/d^j,)^^ (the in- 
dex denotes differentiation at fixed bare parameters). 
From the renormalization of the coupling constant (|^) 
one gets /3 = (— e -I- 87 — 71,) w. At the stable nontriv- 
ial fixed point 7^ with /3 (u*) = (for e > 0), 
this equation yields the anomalous scaling dimension 
T] = ^ (m,) = e/3 -I- 2(5 of the height-field s where I have 
defined 5 = 7ti (w*) /6. From the calculated renormal- 
ization factors (g|) and (|l^), the 7- functions are derived 
as 7 = 2"-2w + O and 7„ = 2Bau'^ + O {u^). Thus 
from [3 (m») = one finds to first order in the e-expansion 

— l?^°'el'i + O (e^). Hence the correction is given by 



27 

_ / 0.01361 (e/2) 
0.002489 £2 



if 
if 



a 
a 



(13) 



For the physical case of a two-dimensional surface with 
nonconserved noise (a = 0, dc = 4, £ = 2) or a one- 
dimensional surface with conserved noise (a = 1, dc = 2, 
£ = 1) the numerical prefactors represent the results 
for 5. As I have stated above, these corrections are al- 
ways small but nonzero positive numbers. From the 
anomalous dimension one finds the roughening exponent 
X = (e — 77) /2 and the dynamic exponent z — 4 — 77 as an- 
ticipated in (H) . Of course the scaling law z — 2% = 4 — e 
holds, being an exact consequence of the trivial renor- 
malization of the noise. But the relation 4 = z -I- x = 
4 -I- (£ — 3?]) /2 = 4 — 35, where the first equality comes 
from the ill-defined symmetry (^), is definitely wrong. 

In summary I have shown that in contrast to a state- 
ment found in most of the recent publications on growth 
models with surface diffusion, the coupling constant of 
these models renormalizes non-trivially. This implies 
that the widely accepted scaling exponents are to be cor- 
rected. But my two-loop result also shows that these 
corrections are very small as long as one accepts the nu- 
merical result at least as a qualitative one with respect to 
the order of magnitude. Thus the supposedly exact expo- 
nents seem to be a very good approximation and the devi- 
ations are not seen in simulations and experiments up to 
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now. I suppose that the origin of the smallness of the cor- 
rection is found in the gradient couphng leading to many 
(but incomplete) cancellations between diagrams as well 
as within the internal momentum integrals. Therefore I 
expect the mode-coupling approach which neglects ver- 
tex corrections completely to be appropriate here. I will 
come back to this in the future. 

I remark that the calculation presented in this letter is 
the first two-loop calculation for the vertex functions of 
the stochastic models considered here. Concerning these 
models the work of Das Sarma and Kotlyar ||2^ con- 
tains a collection of two- loop self-energy diagrams. How- 
ever these diagrams are irrelevant for the problem of the 
second-order renormalization of the coupling constant. 
Thus the authors cite the one-loop result x + ^ = 4 for 
the model with nonconserved noise. Based on the ill- 
defined transformation (||) they believe this scaling law 
to hold to all orders in perturbation theory. As I have 
shown this statement is incorrect. 
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FIG. 1. Cancellation of diagrams: a) an earliest external 
vertex, b) an external vertex inserted in a two-line segment, 
but c) the insertions of an external vertex in a three-line seg- 
ment in all possible ways do not cancel in general, (time in- 
creases always from right to left, the arrows indicate the time 
direction of the propagators). 

FIG. 2. The fundamental diagram: a) an internal earliest 
vertex with all its propagator permutations, b) the two-loop 
kernel diagram, (time and momenta fiow from right to left). 
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Fig. 2a 
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Fig. 2b 
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